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科学 院 行星 科学 重点 


摘要 : 对 上 海天 文 台 火 星 探测 器 定 轨 软 件 的 动力 学 模型 进行 优化 和 升级 ， 开 发 形成 了 MarsODP 
V2.0， 该 软件 可 以 处 理 多 种 类 型 测量 数据 。 对 火星 进行 高 精度 动力 学 建 模 ， 考 虑 火星 本 身 特 
性 ， 增 加 了 火 卫 一 和 火 卫 二 的 三 体 摄 动 、 火 星 大 气 阻力 摄 动 以 及 多 种 火星 定向 模型 优化 等 。 利 
用 MEX 实测 数据 比较 了 优化 前 后 的 火星 定 轨 软件 的 定 轨 精 度 ， 计 算 结 果 表明 ， 动 力学 模型 优化 
， 定 轨 精 度 明显 提升 ， 与 参考 轨道 的 位 置 差异 和 速度 差异 减 小 ， 利 用 实测 数据 对 美国 MRO 探 
测 器 进行 精密 定 轨 并 验证 软件 优化 后 定 轨 精度 。 对 “天 问 一 号 ”探测 器 环 火 轨 道 进行 定 轨 计 算 分 


of 


析 。 
关 键 W: 火星 探测 器 ， 精密 定 轨 ; MRO; MEX; 天 问 一 
中 图 分 类 号 : P173.3 文献 标识 码 : A 

1 引 言 


火星 探测 是 人 类 通过 向 火星 发 射 空间 探测 器 从 而 实现 对 火星 科学 探测 活动 。 探 究 火 星 
及 掌握 火星 的 地 质地 貌 、 水 文大 气 、 资 源 利用 等 方面 ， 都 有 着 极为 重要 的 意义 。 近 现代 ， 人 
类 利用 空间 探测 器 对 火星 进行 深 空 探测 研究 几乎 贯穿 整个 航天 史 。 
美国 在 火星 探测 方面 经 验 丰 富 ， 美 国 发 射 的 火星 全 球 勘 测 者 (Mars Global Surveyor, 
MGS)、 火 星 勘 测 轨道 飞行 器 (Mars Reconnaissance Orbiter, MRO) 等 探测 器 定 轨 精度 可 以 


达到 10 m 左右 ”。 昌 然 中 国 深 空 探测 器 定 轨 工作 在 探 月 工程 中 积累 了 大 量 经 验 “， 但 
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火星 探测 任务 对 中 国 


的 天 文学 家 提出 了 新 的 挑战 。 


火星 的 空间 环境 与 


天 文学 进展 


j 球 也 有 所 不 同 。 如 火星 拥有 两 颗 天 然 


= (Deimos)， 火 星 表 面 存在 稀薄 大 气 "， 火 星 的 重力 场 文件 需 


模型 。 


火星 勘测 轨道 飞行 器 、MAVEN 轨道 器 、 
器 (Mars Orbiter Mission, MOM), 


A 


配 " 等 。 因 此 ， 高 精度 火星 探测 器 精密 


41 卷 


与 月 球 探测 相 比 ， 火 星 距 离 遥 远 ， 并 
卫星 火 卫 一 


(Phobos) 和 火 卫 


要 与 对 应 的 定向 模型 相 匹 


者 密 定 轨 软 件 的 开发 需要 充分 考虑 


足够 精确 的 动力 学 


目前 有 8 颗 探 测 器 在 火星 轨道 中 运行 ， 分 别 是 美国 世 


(T ianwen-1). 火星 


(Curiosity), RAS oe 以 及 中 


火星 表面 上 还 有 


(Ingenuity)。 


已 经 结束 任务 的 旅居 者 号 
号 (Opportunity, MER-B) 和 凤凰 号 (Phoenix); 


2020 年 7 月 


ee 


10 月 9 日 ,“ 天 问 一 


IKERNE RUI” f 


利 完成 深 空 机 动 ， 昨 


ie 


探测 器 于 2 月 24 日 成 功 进入 
该 轨道 上 维 ay 80 d 的 飞 
H7185, KÆ 


5 A174, SKERA i 


(Sojourner), 


星 奥 德 


KE Fa AL 


ak 


i 赛 写 (Odyssey, ODY), 

欧 空 局 的 火星 快车 号 (Mars Express, MEX), EP 
IREA: 
Æ (ExoMars Trace ex Orbiter, TGO)、 阿 联 酋 的 希望 号 (Hope) 以 及 中 国 的 “天 问 一 号 ” 
许多 火星 车 和 着 陆 器 ， 其 中 还 如 
国 的 祝融 号 (Zhurong)， 


勇气 号 


:工作 的 火星 车 有 美国 的 好 帮 号 
和 器 洞察 号 (InSight); 
(Spirit，MER-A)、 机 遇 
除 此 之 外 还 有 首 架 火星 无 人 飞机 机 智 号 


ETE PA SC 


航天 发 射 场 由 长 征 五 号 


中 国 第 一 次 自主 火星 探测 任务 。 经 过 两 次 中 途 修正 ，2020 年 


才 获 ， 成 为 一 颗 环 火 探测 器 ， 
入 近 火 点 高 度 为 280 km 的 大 机 


经 过 远 火 点 


次 经 过 两 次 中 途 人 修正， 于 2021 2 H10 H 
变 轨 和 两 次 近 火 点 制 动 ， 


行 ， 为 火星 车 着 陆 点 区 域 进行 成 
与 环绕 器 分 离 ， 并 成 功 着 陆 于 火星 表面 ， 成 为 中 国 首 颗 火 星 着 陆 器 ; 


维持 火星 车 与 地 球 之 间 的 通 


探测 。 


本 文 首先 在 上 海天 文 台 已 有 火星 探测 器 定 轨 软 件 MarsODP” 基 而 
轨 软 件 的 模型 和 功能 进行 优化 和 升级 ， 


精度 。 


化 后 软件 的 定 轨 精度 进行 比较 分 析 。 
正确 性 和 精度 进行 验证 。 


最 后 对 “天 问 


2 火星 探测 器 动力 学 模型 优化 


针对 我 国 的 自主 火星 探测 任 


一 号 ”探测 器 的 环绕 段 i 


务 “ 天 问 一 bake 


Dea ER 
象 ， 选 取 合适 的 着 陆 点 。5 月 15 


开发 形成 了 MarsODPV2.0, 3 
其 次 通过 对 美国 MRO 探测 器 进行 精密 定 轨 ， 对 软件 的 


入 轨道 周期 为 8.2 h 的 中 继 维持 通信 轨道 (250 km x 12000 km), 
。11 月 8 日 环绕 器 变 轨 至 遥感 使 命 轨 道 ， 开 展 火 星 全 球 遥 感 


上， 对 火星 探测 器 定 


JH MEX 数据 对 优 


轨道 计算 


行 定 轨 计算 ， 并 评估 轨道 


团队 (简称 上 海 台 团 


队 ) 在 已 有 的 火星 探测 器 定 轨 软 件 MarsODP™ 基础 上 ， 对 火星 探测 器 定 轨 软 件 的 模型 和 功 


能 进行 优化 和 升级 ， 开 发 形成 了 MarsODP V2.0。 BAKY 
MarsODP V2.0 定 轨 软件 


在 太阳 系 质心 广义 相对 论 参 考 忆 


所 增 功 能 和 模型 见 表 1。 
匡 架 下 对 观测 量 进 行 高 精度 建 


2 期 单 荃 ， 等 : 火星 探测 器 动力 学 模型 优化 及 “天 问 一 号 ”探测 器 定 轨 研 究 


表 1 


MarsODP 与 MarsODP V2.0 功能 及 模型 比较 


MarsODP 


MarsODP V2.0 新 增 


处 理 的 测量 数据 类 型 


双 程 /三 程 测 距 、 单 程 / 双 程 


测速 、VLBI 时 延 /时 延 率 


光学 测 角 、 同 波束 VLBI 


定 轨 定 位 算法 动力 学 统计 定 轨 运动 学 统计 定位 / 定 轨 
处 理 观 测 目标 环绕 探测 器 火星 着 陆 器 

中 心 天 体 引 力 、 非 球 型 引力 火 卫 一 引力 摄 动 、 火 卫 二 
动力 学 模型 摄 动 、N 体 摄 动 、 光 压 摄 动 、 ”引力 摄 动 、 火 星 大 气 阻 力 


IAU 定向 模型 、 相 对 论 摄 动 


摄 动 、Pathfinder 定向 模型 
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模 ， 在 原 有 处 理 测 距 、 多 普 勒 和 VLBI 时 延 、 时 延 率 测 量 数据 基础 上 ， 增 加 了 光学 测 角 、 同 


波束 VLBI 等 多 种 测量 类 型 ， 


模式 的 计算 功能 。 


具备 动力 学 统计 定 轨 、 运 动 学 统计 
该 软件 除了 处 理 单 目标 绕 飞 探测 器 外 ， 还 可 以 对 火 面目 标 ; 


为 了 提高 定 轨 精度 ， 上 海 台 团 队 进行 了 高 精度 火星 探测 器 动力 学 模型 建 模 ， 帮 有 
力学 模型 基础 上 增加 了 火星 卫星 摄 寻 


2.1 ”火星 探测 器 摄 动力 量 级 估计 
作用 于 卫星 的 某 摄 动力 量 级 通常 采用 摄 动力 与 中 心 引力 之 比 来 估计 ， 可 以 表示 实际 轨 
道 因为 摄 动 力 的 存在 而 偏离 二 体 问 题 轨道 的 程度 。 


设 P 为 摄 动 加 速度 ，X 为 该 摄 动 之 量 级 ， 则 有 


其 中 ，GM 为 中 心 天 体 ( 火 
道 高 度 为 100 ~ 1000 km 的 火星 环绕 卫星 ， 


摄 动 ， 对 于 低 轨 卫星 还 需 考虑 大 气 阻力 的 影响 ， 表 2 给 出 作用 于 火星 探测 器 各 摄 动力 的 量 


级 大 小 * $ 


=I 


E) WAS RCRA ERER PBA. MPH 


«=P (| 


=) 


定 轨 以 及 运动 学 统计 定位 三 
行 定位 。 

E 原 有 
j、 火 星 大 气 阻力 摄 动 以 及 火星 定向 模型 优化 等 。 


的 动 


(1) 


主要 摄 动 源 是 非 球形 引力 ， 其 次 是 第 三 体 引 力 


表 2 作用 于 火星 探测 器 的 各 摄 动 力 之 量 级 大 小 
摄 动 名 称 量 级 
太阳 摄 动 1078 
N 体 摄 动 火星 卫星 摄 动 10~°(Phobos), 107" (Deimos) 
其 他 行星 摄 动 < 10-72 
太阳 直射 光 压 摄 动 < 110? 
J 10-3 
火星 非 球形 引力 摄 动 J2,2 1074 
Cim, Sım, L23 1075 ~ 107° 
火星 大 气 阻力 摄 动 万 二 200 a 
相对 论 摄 动 10-10 
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综 上 分 析 ， 本 文 定 轨 分 析 时 使 用 上 海天 文 
力学 模型 及 定 轨 部 分 策略 见 表 3。 考 虑 了 火星 
大 气 摄 动 。 其 中 NN 体 摄 动 除了 考虑 了 太 民 
卫星 的 火星 而 言 ， 需 要 考虑 火 ] 


FP 心 引力 ， 


{TA 


ERA 


火星 探测 器 定 轨 动 力学 模型 


月 引力 摄 动 外 ， 对 于 
D 二 的 三 体 摄 动 以 提高 定 轨 精 度 。 


41 卷 


火星 定 轨 软 件 ， 对 火星 探测 器 定 轨 考 虑 的 动 
非 球 型 引力 摄 动 ，N 体 摄 动 、 火 星 


TATAR 


项 目 Li 
基本 参考 系 KORRES AR 
中 心 天 体 火星 
非 球形 引力 摄 动 重力 场 模型 MRO120D; 采用 pathfinder 定向 模型 
N 体 摄 动 太阳 、 大 行星 、 火 卫 一 、 火 卫 二 (采用 
太阳 辐射 压 固定 面 质 比 模型 
火星 大 气 摄 动 指数 模型 
数据 使 用 及 权重 设置 ” 测 距 为 1 m， 测 速 为 1 mm/s， 时 延 为 0.1 ns， 时 延 率 为 1 ps/s 
大 气 延 迟 修正 Saastamoinen-NMF 
a} at KSG 二 阶 定 步 长 积分 器 


2.2 ”火星 探测 器 动力 学 模型 精 化 


相 较 于 月 球 探测 ， 火 星 探测 需要 考虑 火 司 
定向 模型 选取 的 影响 。 本 节 主 要 从 以 下 三 个 方面 对 软件 模型 优化 进行 曾 述 。 
221 火 卫 摄 动力 模型 建立 


此 前 上 海天 文 台 综合 定 轨 软件 
而 对 火星 探测 器 的 精密 定 轨 考 虑 火电 

火 卫 一 (Phobos) 与 火 卫 二 (Deimos) 
以 近 圆 形 的 轨道 公转 。 火 卫 一 是 两 颗 卫 星 中 更 大 的 一 颗 ， 三 轴 半 径 为 (13.3 x 11.1 x 9.3) + 
0.3 km， 质 量 为 10.6x 101 kg， 密 度 为 1900 kg/m; 火 卫 二 的 三 和 
0.5 km， 质 量 为 2.4 x 1015 kg， 密 度 为 1750 kg/m3。 火 星 ] 
火 卫 一 绕 火 星 运动 半 长 轴 为 9376 km， 围 绕 火星 一 周 
行 的 距离 为 23 459 km， 轨 道 周 期 为 1.26244 d' 


PN 体 摄 动 只 考虑 了 太 B 


:卫星 的 第 三 


卫星 摄 动 、 火 星 大 气 阻力 摄 动 以 及 不 同 的 火星 


系 中 的 大 行星 和 
体 引力 


区 状 不 规则 并 充 ; 


pitt 


A A 2 ARIE A Ak 


AS SCT AFD Ak FP A 


有 星 历 表 为 JPL 在 20 


1877—2011 年 期 间 所 有 基于 地 球 和 航天 器 的 观测 ， 包 括 MEX ZRK 
跟踪 数据 和 2011 年 之 前 对 火 了 


Vii 


半径 为 (7.6x 6.2 x 5.4) + 


ae 
只 需要 


二 的 额外 成 像 观测 数据 ， 是 


B 星 虽然 很 小 ， 但 十 分 接近 火星 ， 
0.31891 d; KI 


月 引力 摄 动 ， 


让 


Hl SEK BS 


14 年 发 布 的 MAR097" , REST 
各 火 卫 一 的 多 普 勒 


H 


目前 


MEX 探测 器 的 半 长 轴 约 为 8700 km (火星 半径 3397 km), 4 


星星 历 


高 心率 e 约 为 0.57， 是 一 颗 


大 偏心 率 轨道 卫星 ， 其 近 火 点 高 度 约 300 km、 远 火 点 高 度 超过 1 x 104 km， 轨 道 周期 约 为 


6.9 h。 选 取 MEX PRIA WH 


卫 一 最 近 距 离 为 4538 km. 
影响 比较 大 "。 


[的 时 间 为 2009-08-07T20-00-00 (UTC)， 一 个 了 


于 MEX 为 大 椭圆 


卫 一 和 火 卫 二 距离 火星 的 高 度 如 图 1 所 示 ，MEX 探测 器 到 两 颗 卫 / 
轨道 ， 会 飞 掠 火 卫 


PEHIA, K 
= 2, 距离 火 
星 引力 摄 动 


2 期 HE, SE: 火星 探测 器 动力 学 模型 优化 及 “天 问 一 号 ”探测 器 定 轨 研 究 


距离 /km 


YE: 时 间 自 2009-08-07T20-00-00 起 。a) 火 卫 一 ; b) KEZ. 


图 1 火 卫 一 和 火 卫 二 到 火星 表面 的 距离 变化 示意 图 


4 
1.8*10 
1.4 
g 
wa 
ae 
= 10 
0.6 
20 22 24 
t/h 
a) 


4 
2.820 


2.4 


26 28 20 22 


YE: 时 间 自 2009-08-07T20-00-00 起 。a) 火 卫 一 ; b) KEZ. 


24 26 28 


图 2 MEX 探测 器 到 火 卫 一 和 火 卫 二 距离 变化 示意 图 


2.2.2 火星 定向 模型 优化 


火 固 系 的 定义 需要 与 定 轨 使 用 的 重力 场 模型 相 匹 配 ， 目 前 火星 探测 器 定 轨 软 件 
的 是 TAU 火星 定向 模型 ， 定 轨 使 用 的 是 对 应 于 pathfinder 定向 模型 的 如 
MGS95J 和 MRO120D 等 ， 模 型 不 匹配 会 对 轨道 确定 和 预报 带 来 一 定 的 误差 。 火 星 固 
标 系 与 火星 天 球 参考 系 转换 关系 参照 图 3。 


EE 力 场 模型 
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中 使 
d 如 
BEA 


火星 重力 场 模型 可 以 在 行星 数据 系统 (Planetary Data System, PDS) 网 站 https://pds. 
nasa.gov/ 中 下 载 ， 重 力 场 模型 对 应 的 说 明文 档 介 绍 了 Pathfinder 定向 模型 所 需要 的 参数 。 

MGS95J 中 本 初子 午 面 选取 与 IAU 2000 标准 的 J2000.0 历 元 的 本 初子 午 面 一 致 "， 将 
MGS95 重力 场 模型 对 应 的 Pathfinder 定向 参数 模型 与 IAU 2009 定向 模型 进行 比较 ， 分 别 


使 用 两 个 模型 计算 自转 轴 > 轴 


(0 km, 0 km, 3396 km) 与 z # (3396 km, 0 km, 0 km) 转换 


至 火星 天 球 参 考 系 的 坐标 ， 


首 比 较 两 者 的 差异 。 图 4 给 出 了 两 个 名 


转换 坐标 后 的 差异 : 


从 
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火星 本 初子 午 


火星 天 球 参 考 系 
(ICRF) 


图 3 


2010 一 2030 Æ, x 轴 


ME Ho 


z 轴 差异 /m 


E: 时 间 自 2010 年 起 。 


ay 


图 4 


2.2.3 火星 大 气 密度 模型 建立 
对 比 月 ] 
(95.32%), 


=| 
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x 轴 差 异 /m 
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火星 历 元 真 赤道 
火星 历 元 平 赤 道 


面 


ICRF 基 本 平面 
近似 为 J2000.0 地 球 平 赤道 


X=AB, w=AC, ¢=CD, Q=yB 


火星 固 联 坐 标 系 与 火星 天 球 参 考 系 转换 关系 示意 图 


转换 差异 近 100 m, z 轴 转 换 差 异 近 50 m， 且 周年 与 半 周 年 变化 项 
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不 同 火星 定向 模型 坐标 转换 差异 比较 


球 表面 的 真空 状态 ， 火 星 上 有 稀薄 的 大 气 ， 火 星 大 气 层 的 主要 成 分 是 CO 
次 是 N、Ar， 此 外 还 有 少量 的 O 和 水 燕 气 ， 平 均 压 力 仅 为 700 Pa。 火 星 大 气 


分 为 : 外 逸 层 、 上 层 大 气 (A 


R) 


上 的 部 分 ， 大 气 边界 不 明显 ， 


由 于 火星 大 气 的 存在 ， 对 低 轨 火 星 卫 星 以 及 并 


中 层 大 气 和 低层 大 气 ， 其 中 外 逸 层 通常 是 指 200 km 及 以 
它 只 是 逐渐 消失 ”。 


陆 器 着 陆 阶段 受到 的 大 气 阻力 摄 动 不 可 


忽略 ， 大 气 阻力 摄 动 加 速度 是 关于 火星 大 气 密度 的 函数 ， 火 星 大 气 密度 可 以 通过 火星 大 气 数 


“faj 


据 库 或 火星 大 


化 模型 获得 。 KA 


NAR T KEKEE, KAM S KARE ESA, 


耗 时 过 长 ， 不 适用 于 实时 仿真 。 


火星 大 气 简 化 模型 则 是 利 月 


H 


大 气 数据 库 可 以 


AE 


Hh BQ EH kK EK URE, 8 
庞大 的 数据 资料 使 得 调用 数据 库 


昌 线 拟 合 方式 ， 将 实际 任务 数据 或 火星 大 气 数据 库 提 供 的 


数据 简化 为 指数 模型 ， 其 大 气 信息 比较 单一 ， 模 型 多 为 大 气 密度 随 海拔 高 度 变化 的 一 维 指数 


2 期 


模型 或 分 层 指 数 模 型 


气 数据 库 提供 
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’ Lh ke ifn. 


压强 
下 JERS 


密度 随 高 度 变化 


的 模型 ， 由 实际 任务 数据 或 火星 大 


的 数据 拟 合 得 到 。 模 型 多 为 大 气 密度 随 火 四 


指数 模型 ， 以 及 


中 拥有 计算 量 少 、 


火星 大 气 模型 ， 


Ulin. 压 强 、 
速度 快 的 优点 ， 


Ph = po , po =1.58x 107? kg/m? , h, = 9354.5m 


Joel Benito 2008 一 维 指数 模型 
器 的 位 置 即 可 求 出 当前 的 密度 ， 


密度 随 高 度 变 化 的 模型 ” 
适合 在 实时 定 轨 和 仿真 分 析 中 使 用 。 
选择 的 是 简单 的 一 维 指数 模型 *”， 其 
在 实时 定 轨 和 仿真 分 析 中 使 用 。 


高 度 变化 的 一 维 指数 模型 或 分 层 
。 火 星 大 气 简化 模型 在 实际 应 用 
式 (2) 为 本 文 使 用 的 
速度 快 的 优点 ， 适 合 


20. 


\ 有 计算 量 少 、 


(2) 


» po 火星 表面 的 大 气 密 度 ，j。 均 质 大 气 高 度 ，h 为 探测 器 的 高 度 ，pn 为 所 求 高 度 的 大 


2.3 ”软件 定 轨 精 度 比 较 


2009 年 8 


VLBI 站 以 及 欧 空 局 新 诺 舍 深 空 站 组 织 了 对 MEX 的 观测 试验 。 新 诺 舍 深 空 站 在 X 频段 发 射 


17 日 20:00 至 次 


上 行 信号 ， 经 星 上 转发 器 转发 后 ， 分 为 8 和 XX 两 个 频段 的 下 行 信号 ， 由 国内 测 站 接收 其 下 
行 信号 ， 分 别 在 S 频段 和 X 频段 上 获取 了 三 程 多 普 勒 数据 (5 s 积分 时 长 )， 同 时 还 获得 
空 局 NNO 测 站 跟踪 MEX 的 双 程 多 普 勒 数据 (1 s 积分 时 长)。 

利用 该 弧 段 的 MEX 数据 ， 分 别 选择 不 同 的 软件 版 本 进行 定 轨 计 算 分 析 ， 定 轨 结 果 与 


日 04:13 (UTC)， 


中 大 气 密度 只 与 探测 器 到 火星 表面 距离 有 关 ， 


再 通过 公式 求 得 探测 器 受到 的 火星 大 气 阻力 加 速度 。 


民 据 探测 


国内 测 站 上 海 、 昆 明 、 马 鲁 木 齐 3 个 


TEK 


Kk 


空 局 提供 的 参考 轨道 (比利时 皇家 天 文 台 精 密 轨道 ， 精 度 约 20 ~ 25 m”) 进行 星 历 比较 ， 
轨道 差异 见 表 4 
表 4 火星 卫星 摄 动 对 MEX 定 轨 精 度 影响 
ee (eS ae 
MarsODP 1.78 357.64 162.02 392.64] 0.06 0.029 0.049 0.082 | 0.14 0.10 
MarsODP V2.0 | 1.04 6.99 3.80 8.02 | 0.001 0.001 0.001 0.001 | 0.13 0.10 


计算 结果 


米 下 降 到 10 m 左右 ， 速 度 差 异 从 约 0.1 m/s | 


明 ， 优 化 后 的 软 伯 


F 定 轨 精 度 有 显著 提高 ， 与 参考 轨道 的 位 置 差异 从 300 多 


下 降 至 0.001 m/s。 差 异 较为 明显 主要 原因 


KE 


MEX 探测 器 为 大 椭圆 轨道 ， 与 火星 卫星 距离 较 近 ， 火 星 卫 星 摄 动 对 MEX 速度 影响 明显 ， 


=] 


在 加 入 火星 


卫星 


摄 动 后， 测速 观 测量 可 以 更 好 二 


3 MRO 探测 器 实测 数据 定 轨 分 析 


美国 航天 局 


4124 


RHIZ, 


也 约束 MEX 轨道 ， 定 轨 精 度 明显 提升 。 


的 火星 勘测 轨道 飞行 器 (Mars Reconnaissance Orbiter, MRO) 于 2005 年 8 


并 于 2006 年 3 


进入 火星 轨道 。 
和 地 表 条 件 ， 探 测 地 表 下 的 水 冰 和 矿物 ， 监 测 日 常 天 气 ， 并 为 未 来 的 任务 确定 着 


MRO 配备 了 一 套 先进 的 仪器 ， 用 于 看 
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陆地 点 。 轨 道 器 还 测试 了 一 个 新 的 通信 系统 ， 该 系统 能 够 以 更 快 的 速度 向 航天 器 传输 数据 ， 
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使 MRO 能 够 作为 其 他 任务 的 重要 中 继 卫 星 。MRO 飞行 轨道 高 度 为 255 km x 320 km 圆 轨 


道 ， 轨 道 周期 约 2 h”。 


本 节 对 MRO 探测 进行 定 轨 精 度 分 析 ， 观 测 数据 为 美国 DSN 深 空 站 双 程 多 普 勒 测量 数 
据 。DSN 深 空 站 在 全 球 布 网 且 分 布 均匀 ，2017 年 1 月 至 4 月 期 间 多 个 DSN 深 空 站 对 MRO 
接续 跟踪 观测 ， 每 天 观测 12 ~ 24 np， 观测 弧 段 分 布 均匀 。 


利用 2017 年 1 月 至 4 月 的 双 程 多 普 勒 测量 数据 进行 定 轨 分 析 ， 定 轨 使 用 的 动力 学 模型 
参数 设置 见 表 3， 定 轨 弧 长 24 h， 积 分 步 长 10 s， 同 时 解 算 光 压 系数 和 大 气 阻 力 系数 。 四 个 


弧 段 的 定 轨 结 果 与 NASA 发 布 的 重建 轨道 进行 比较 。 由 图 5 可 以 看 出 1 月 份 解 算 的 轨道 差 
异 较 大 ， 之 后 轨道 的 差异 逐渐 变 小 ， 位 置 差异 小 于 30 m、 速 度 的 差异 小 于 0.03 m/s。 


位 置 差异 /m 


Pos/m 


RI 


T 方 向 


速度 差异 /(m-s!) 


入 方向 


V.,/(ms") 


TE: 时 间 从 2007-01-01 起 。 


5 2017 年 1 一 4 月 解 算 轨 道 与 NASA MRO 重建 轨道 位 置 与 速度 差异 


图 6 给 出 2017 年 1 月 至 4 月 不 同 定 轨 精 度 所 对 应 的 天 数 之 和 : 与 NASA 的 MRO 重建 


轨道 位 置 差异 小 于 10m HFA 36d. 10m 至 20 m WA 38d, 20 m 至 30 m 的 有 23 d; 速 


度 差 异 小 于 0.01 m/s 的 40 d、0.01 m/s 至 0.02 m/s 的 有 39 d、0.02 m/s 至 0.03 m/s 的 有 


20 do 
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40 


30 50 70 90 j 0.01 0.03 0.05 0.07 0.09 


位 置 差异 /m 速度 差异 / (ms!) 
图 6 2017 年 1 一 4 月 解 算 轨 道 与 NASA MRO 重建 轨道 差异 统计 图 


从 2017 年 1 一 4 月 数据 定 轨 结 果 中 ， 每 月 各 选取 1 d 弧 段 定 轨 结果 ， 其 星 历 比较 结果 见 
表 5， 计 算 结 果 表 明 , 位 置 误 差 主要 在 丰 和 方向 上 ，RR 方 向 误差 小 于 0.7 m， 测 速 数据 残 
差 在 0.6 mm/s 左右 。 


表 5 MRO 探测 器 定 轨 星 历 比较 结果 


定 轨 日 其 位 置 /m 速度 /(cm.s-!) ae 
R T N Pos R T N Ve | /(mm:s ) 
2017-01-26 | 0.68 27.14 12.03 29.69 | 2.53 0.06 1.12 2.77 0.59 
2017-02-23 | 0.59 4.15 2.02 4.65 | 0.38 0.05 0.19 0.42 0.50 
2017-03-23 | 0.51 10.76 6.31 12.49 | 1.00 0.05 0.59 1.16 0.52 
2017-04-05 | 0.59 2.68 1.03 2.93 | 0.23 0.06 0.10 0.26 0.66 


综 上 所 述 ， 使 用 2017 年 1 一 4 月 的 单 日 双 程 多 普 勒 数据 定 轨 ，1 s 积分 时 长 的 斜坡 双 程 
多 普 勒 数据 定 轨 后 残 差 约 为 0.5 ~ 0.6 mm/s， 解 算出 的 轨道 与 NASA 重建 轨道 位 置 差异 小 
于 30 m， 速 度 的 差异 小 于 0.03 m/s. 


4 “天 问 一 号 ” 环 火 轨道 定 轨 精 度 分 析 


本 章 主 要 对 火星 环绕 器 进行 精密 定 轨 ， 通 过 重 有 弧 段 比较 方法 分 析 轨 道 确 定 精度 。 
4.1 “天 问 一 号 ”停泊 轨道 

2021 年 2 月 24 日 第 三 次 近 火 制 动 后 ， 火 星 探测 器 进入 50 h 周期 的 大 椭圆 停泊 轨道 ， 
近 火 点 高 度 为 280 km， 远 火 点 高 度 约 为 59 000 km。 环 绕 器 通过 星 上 多 种 科学 载荷 对 火星 空 
间 大 气 、 气 象 、 地 形 等 进行 科学 探测 ， 绘 制 高 精 度 着 陆 区 域 表 面 的 地 形 图 ， 选 取 合适 的 着 陆 
点 。 为 了 保证 科学 研究 的 顺利 展开 ， 尤 其 是 对 火星 着 陆 区 域 的 高 精度 成 像 ， 需 要 高 精度 的 探 
测 器 轨道 信息 。 在 停泊 轨道 飞行 过 程 中 ， 三 个 深 空 站 每 天 接力 进行 观测 ，VLBI 观测 网 每 两 
天 进行 一 次 观测 ， 以 满足 高 精度 的 测定 轨 需 求 。 
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利用 2021 年 2 月 24 日 一 5 月 15 日 之 间 的 观测 数据 进行 ; 
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情 密 定 轨 ， 定 轨 弧 长 约 3 ~ 5 


圈 (6 ~ 10d), EIK 1 圈 (2 d)， 测 量 数据 使 用 测 距 测速 联合 VLBI 数据 ， 重 登 精度 分 


MERWE 7 所 示 。 定 轨 后 残 差 ， 测 距 约 为 0.5 m, ME 


为 0.16 mm/s， 时 延 为 0.1 ns, H} 


延 率 为 0.33 ps/s。 


302 318 


300 2 


位 置 误差 /m 


6.0 54 55 


4.0 


2.0 


速度 误差 /(m.s7) 


图 7 


计算 结果 表明 ， 对 于 火星 停泊 
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Heme 


利用 联合 


轨道 ， 


火星 停泊 轨道 探测 器 重印 弧 段 定 轨 分 析 位 置 速度 误差 


距 测 速 和 VLBI 数据 进行 计算 可 得 ，10 


= 


WUE 


d 弧 长 定 轨 位 置 精度 好 于 300m, HAR 方向 好 于 120 m， 误 差 主要 在 人 方向 上 ， 速度 精 


度 好 于 0.035 m/s. 
4.2 “天 问 一 号 ”维持 轨道 
2021 年 5 H 


15 日 7 时 18 分 火星 车 成 功 着 陆 


火星 乌托邦 平原 ，49 h 后 环绕 器 进行 轨 


控 ， 进 入 轨道 周期 为 8.2 h 的 中 继 通信 维持 轨道 ， 近 火 点 高 度 为 250 km， 远 火 点 高 度 约 为 


12000 km. 
可 实现 一 个 火星 日 
命 轨 道 ， 开 展 火星 全 球 遥 感 探 测 。 


中 继 轨 道 为 每 个 火星 日 
内 对 火星 车 的 近 火 和 远 火 两 次 中 继 通信 。 


运行 3 


5 月 30 日 之 前 为 密集 观测 ， 三 个 深 空 站 对 环绕 
5 月 30 日 之 后 为 常规 观测 ， 由 双 站 进行 接力 观测 ， 


图 的 回归 轨道 ， 每 个 火星 日 都 会 经 过 着 陆 点 上 空 ， 


11 月 8 日 环绕 器 变 轨 至 遥感 使 


VLBI 每 天 观测 4 h; 
一 周 观 测 三 


器 进行 接力 观测 ， 
每 天 观测 弧 段 缩短 ，VLBI 


次 ， 每 次 1 h。 每 次 轨道 机 动 前 后 会 进行 一 段 时 间 的 密集 观测 。 
本 节选 取 了 2021 Æ 5 H 24 HZ6 H11 日 “天 问 一 
度 分 析 。 该 弧 段 轨道 小 机 动 频繁 (8 ~ 12h 
距 联 合 VLBI 数据 ， 并 分 弧 段 解 算 测 距 系 统 误差 。 
HAI 48 h 左右 弧 段 (5 ~ 7 圈 ) HA 1 (8h), 


号 ”环绕 器 的 实测 数据 进行 定 轨 精 
次 ) ， 而 需要 解 算 多 组 经 验 力 。 测量 数据 使 用 测 


如 图 8 所 示 ， 定 轨 后 测 距 数 据 
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残 差 好 于 0.4 m，VLBI 时 延 数据 残 差 好 于 0.2 ns. 


05-25T00:00 05-25T12:00 05-26T00:00 
时 刻 (UTC, 2021 年 ) 


时 延 残 差 /ns 


05-25T12:00 05-25T22:00 05-26T08:00 05-26T18:00 
时 刻 (UTC, 2021 年 ) 


8 星 维持 轨道 弧 探测 器 测量 数据 定 轨 后 残 差 图 


轨道 重 又 弧 段 分 析 比 较 结果 见 图 9， 计算 结果 表明 ， 维 持 轨 道 段 定 轨 平均 位 置 精度 好 于 
180 m (其 中 径 向 好 于 3m), RAEBAT AN 方向 上 ， 速 度 精度 好 于 0.035 m/s; 密集 观 
测 弧 段 定 轨 精度 较 高 ， 位 置 精 度 平均 好 于 100 m， 常 规 观测 弧 段 定 轨 精度 较 差 ， 平 均 位 置 精 
度 约 为 350 mo 
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本 文 针 对 我 国 火星 探测 任务 ， 完 善 并 优化 了 已 有 的 上 海天 文 台 火 星 探 测定 轨 软 件 ， 增 加 


了 火星 卫星 摄 动 、 火 星 定 向 模型 优化 以 及 火星 大 气 摄 动 。 


利用 优化 后 的 高 精度 火星 定 轨 软 件 对 MRO 探测 器 进行 轨道 确定 ， 与 参考 轨道 比较 ， 


轨道 平均 位 置 精度 优 于 30 m， 速 度 的 精度 优 于 0.03 m/s， 验 证 了 定 轨 软件 的 正确 性 与 可 


于 300 m， 其 中 有 R 方 同好 于 120 m， 误 差 主要 在 人 方向 上 ; 速度 精度 好 于 0.035 m/s。 维 持 
轨道 定 轨 平 均 位 置 精度 好 于 180 m (其 中 径 向 好 于 3 m)， 速 度 精 度 好 于 0.035 m/s; 当 密 集 


山 
m 
Ly 
oF 


对 “天 问 一 号 ”探测 器 的 火星 停泊 轨道 进行 轨道 确定 ，3 ~ 5 圈 弧 长 定 轨 位 置 的 精度 好 


观测 时 ， 定 轨 精 度 较 高 ， 位 置 精 度 平均 好 于 100 m 
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Dynamic Model Optimization of Mars Probe and Orbit 


Determination Analysis of Tianwen-1 


SHAN Quan! ”, LI Pei-jial34, FAN Min®, HUANG Yong!?34+, YANG Peng!? 
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Abstract: In this paper, the dynamic model of Mars probe orbit determination software 
of Shanghai Observatory is optimized and upgraded, which can process various types of 
measurement data and carry out high-precision dynamic modeling of Mars. Considering 
the characteristics of Mars, the three body perturbation model of Phobos and Demos, the 
Mars atmospheric drag perturbation and the optimization of various Mars pointing models 
are added. The MEX orbit determination results using different softwares are compared. 
The calculation results show that the orbit determination accuracy is significantly improved 
after the optimization of dynamic models, the position difference decreased from more than 
300 m to about 10 m, and the velocity difference decreases from about 0.1 m/s to 0.001 m/s; 
The orbit determination result of MRO probe shows that the residual error of the slope 
two-way Doppler data with 1 s integration time is about 0.5 ~ 0.6 mm/s, and the positon 
difference between the calculated orbit and the NASA reconstructed orbit is less than 30 m, 
and the difference of velocity is less than 0.03 m/s; The orbit determination of parking orbit 
of Tianwen-1 probe shows that the position accuracy of 6 ~ 10 d arc length can be better 
than 300 m, in which the R direction is better than 120 m, and the error is mainly in the 
T direction; The velocity accuracy is better than 0.035 m/s, which can meet the scientific 


research index. 


Key words: Mars probe; orbit determination; MRO; MEX; Tianwen-1 


